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ABSTRACT 

We study the effects of initial and boundary conditions, taking two-dimensional hy- 
drodynamical numerical simulations of hot accretion flow as an example. The initial 
conditions considered include a rotating torus, a solution expanded from the one- 
dimensional global solution of hot accretion flows, injected gas with various angular 
momentum, and the gas from a large-scale numerical simulation. Special attention is 
paid to the radial profiles of the mass accretion rate and density. Both can be described 
by a power-law function, M oc r s and p oc r~ p . We find that if the angular momentum 
is not very low, the value of s is not sensitive to the initial condition and lies within a 
narrow range, 0.47 < s < 0.55. However, the value of p is more sensitive to the initial 
condition and lies in the range 0.48 < p < 0.8. The diversity of the density profile 
is because different initial conditions give different radial profiles of radial velocity 
due to the different angular momentum of the initial conditions. When the angular 
momentum of the accretion flow is very low, the inflow rate is constant with radius. 
Taking the torus model as an example, we have also investigated the effects of inner 
and outer boundary conditions by considering the widely adopted "outflow" boundary 
condition and the "mass flux conservation" condition. We find that the results are not 
sensitive to these two boundary conditions. 

Key words: accretion, accretion discs - hydrodynamics- black hole physics 
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1 INTRODUCTION 

Hot accretion flows, such as advection-dominated accretion 
flows (ADAFs; Ichimaru 1977; Rees et al. 1982; Narayan & 
Yi 1994; 1995; Abramowicz et al. 1995; see Narayan, Ma- 
hadevan & Quataert 1998 and Kato, Fukue & Mineshige 
1998 for reviews) are interesting because they are likely op- 
crating in low- luminosity active galactic nuclei (AGNs), in- 
cluding the supermassive black hole in our Galactic center, 
and hard and quiescent states of black hole X-ray binaries 
(see Narayan 2005; Yuan 2007; Narayan & McClintock 2008; 
Ho 2008; Yuan 2011 for reviews). Because of its widespread 
applications, the properties of hot accretion flows have been 
investigated intensively in the last decade. Hot accretion 
flows are originally proposed and studied analytically by 
vertically- integrated one-dimensional method (e.g., Narayan 
1994). In the analytical works, the mass accretion rate is 
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usually assumed to be constant with radius. Under this as- 
sumption, the radial profile of density satisfies p(r) oc r~ 3 ^ 2 . 

Later on, intensive numerical simulation works have 
been performed to investigate the hot accretion flows. One 
of the most important finding is that the mass inflow 
rate is not a constant of radius; rather, it decreases in- 
ward. Correspondingly, the density profile flattens. The den- 
sity profile obtained by the simulations is p oc r~ p , with 
0-5 < p < 1- This is true for both hydrodynamic (HD) 
and magnetohydrodynamic (MHD) simulations (e.g., Stone, 
Pringle & Begelman 1999, hereafter SPB99; Igumenshchev 
& Abramowicz 1999, 2000, hereafter IA99, IA00; Stone & 
Pringle 2001; Hawley, Balbus & Stone 2001; Machida, Mat- 
sumoto & Mineshige 2001; Hawley & Balbus 2002; Igu- 
menshchev, Narayan & Abramowicz 2003; Pen, Matzener & 
Wong 2003; De Villiers, Hawley & Krolik 2003; Yuan & Bu 
2010; Pang, et al. 2011; McKinney, Tchekhovskoy & Bland- 
ford 2012; Li, Ostriker & Sunyaev 2013; see Yuan, Wu & 
Bu 2012 for a review). One model proposed to explain the 
above result is ADIOS (adiabatic inflow-outflow solution). 
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In this scenario, the inward decrease of inflow rate is because 
of mass lost via outflow (Blandford & Begelman 1999, 2004; 
Begelman 2012). 

By analyzing the HD and MHD numerical simulation 
data, Yuan, Bu & Wu (2012, hereafter YBW2012; see also 
Li, Ostriker & Sunyaev 2013) show that outflow must exist, 
thus confirming the ADIOS scenario. Moreover, they investi- 
gated the physical mechanism of producing outflows. In the 
HD case, they find that the temperature of outflow is sys- 
tematically higher than that of inflow (Fig. 8 in YBW2012). 
This suggests that the outflow is driven by the buoyancy 
which arises because of the convective instability of the HD 
accretion flow. In the case of MHD flow, they find that the 
specific angular momentum of outflow is close to the Keple- 
rian angular momentum; while that of inflow is much lower 
(Fig. 11 in YBW2012). This suggests that it is the centrifu- 
gal force associated with the magnetic field that drives the 
outflow (Blandford & Payne 1982). 

The radial profile of mass accretion rate (or density) is 
important because of the following reasons. The profile of 
mass accretion rate determines the emitted spectrum and 
other radiative features of an accretion flow (e.g., Quataert 
& Narayan 1999; Yuan, Quataert & Narayan 2003). With 
the increasing power of telescopes, the gas properties (den- 
sity, temperature) at the Bondi radius of galaxies can been 
observed directly in some sources, such as Sgr A* and 
NGC 3115. Using the observational data and Bondi accre- 
tion theory (Bondi 1952), we can obtain the mass accretion 
rate at Bondi radius. Then depending on the radial profile 
of mass accretion rate (or more exactly the density), the 
radiation of the accretion flow is completely different. In 
cosmology simulations and other simulations studying the 
accretion flow on relatively large scale (e.g., Cuadra et al. 
2005; 2006; 2008; Springel et al. 2005; Booth & Schaye 2009; 
Power et al. 2011; Hobbs et al. 2012), we can at most re- 
solve the Bondi radius and determine the Bondi accretion 
rate. Then the mass accretion rate profile can determine the 
evolution of black hole mass and spin. The radial profile of 
mass accretion rate is also important to determine the feed- 
back efficiency of an AGN. It is well known that feedback 
from AGN can strongly influences the physical properties of 
their vicinity, their host galaxies, and even of the intergalac- 
tic material of galaxy clusters to which they belong (e.g., 
Proga, Stone & Kallman 2000; Proga 2007; Kurosawa & 
Proga 2009; Ciotti & Ostriker 1997; 2001; 2007; Ciotti, Os- 
triker & Proga 2009; Novak, Ostriker & Ciotti 2011). Both 
the radiation and mass outflow from a AGN can influence its 
surroundings. The radiation from AGN can heat or cool the 
gas surrounding it, and the efficiency of heating or cooling 
depends on the spectrum. The outflows can directly blow 
away the gas surrounding an AGN. Thus, if the accretion 
rate at the Bondi radius of a galaxy can be obtained directly 
from observations, the mass accretion rate (or density) pro- 
file can exactly determine how much radiation and outflow 
can be produced and further influence the environment of 
AGN. 

In this paper, we will study how the properties (es- 
pecially radial profiles of mass accretion rate and density) 
of hot accretion flows depend on the initial and boundary 
conditions. Our motivations are as follows. In the realistic 
universe, the environment of accretion flows is likely very 
complex and differ from each other. Correspondingly, the 



initial conditions of accretion flows, such as their radial ve- 
locity, temperature, and angular momentum, should be very 
different in different sources. In the literatures, torus and 
injection models are usually adopted as the initial condi- 
tion of simulations. But no body has studied whether and 
how the different initial conditions affect the properties of 
the accretion flow. YBW2012 studied several HD models. In 
their model A, the initial condition is a rotating torus. In 
model C, they inject gas at the outer boundary. Although 
the description and strength of viscosity are same in differ- 
ent models, the properties of the accretion flows are found 
to be quite different. First, in model A, the flow is nearly 
symmetric to the equatorial plane, with inflow being close 
to the equatorial plane (8 — 70° — 90° in a spherical coordi- 
nate) and outflow close to the surface of the disk (0 = 50° 
and 140°). In model C, the flow is not symmetric to the 
equatorial plane, and large-scale bulk motions dominate the 
flow. Second, in model C, the radial velocity of inflow close 
to the equatorial plane can be described as v r tx r~ ' 5 ; but 
in model A, the radial velocity scaling with radius is much 
steeper. The difference between the results obtained by dif- 
ferent initial settings motivates us to systematically study 
the influence of the initial conditions on the properties of 
hot accretion flows. 

In the present work, we focus on HD simulation. In re- 
ality, people now believe that it is the MHD turbulence as- 
sociated with the magneto-rotational instability (MRI; Bal- 
bus & Hawley 1991; 1998) that is responsible for the angu- 
lar momentum transfer in accretion flows. Therefore, MHD 
simulation is more realistic. Our considerations of adopt- 
ing HD simulation are as follows. First, MHD simulation 
is obviously more expensive; and the results depend on the 
initial configuration and strength of the magnetic field. Sec- 
ond, in many astrophysical simulations such as large-scale 
cosmological ones, magnetic field is not considered and HD 
simulations are still the dominant ones. It is thus of practi- 
cal usefulness to consider HD simulations. Third, we expect 
that if the initial and boundary conditions are important to 
determine some properties of HD accretion flow, it should 
also be important for MHD flows. Moreover, by comparing 
all relevant HD and MHD numerical simulations, Yuan, Wu 
& Bu (2012) found that the radial profiles of inflow rate 
obtained by various HD and MHD simulations are very sim- 
ilar. In other words, it seems that the result does not depend 
on the presence or absence of magnetic field, the strength 
and configuration of the magnetic field, and the dimension 
of calculation (two or three dimension). Such an apparently 
surprising result was predicted and explained in Begelman 
(2012). This implies that maybe the results obtained in the 
present study can be applied to MHD case as well. 

The structure of the paper is as follows. In §2, we de- 
scribe the basic equations, various initial and boundary con- 
ditions adopted and our numerical method. The results of 
simulations will be given in §3. We summarize our results in 
§4. 

2 NUMERICAL METHOD 
2.1 Equations 

In spherical coordinates (r,9,4>), we solve the following hy- 
drodynamical equations describing accretion: 
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dp 
~dt 

dv 



+ pV ■ v = 0, 

-Vp - pVtp + V ■ T, 



r dt 



T 2 /p, 



(1) 
(2) 
(3) 



Here, p, p, v, ip, e and T are density, pressure, velocity , 
gravitational potential, internal energy and anomalous stress 
tensor, respectively. d/dt(= d/dt + v ■ V) denotes the La- 
grangian time derivative. We adopt an equation of state of 
ideal gas p = (7 — l)e, and set 7 = 5/3. 

We use the stress tensor T to mimic the shear stress, 
which is in reality magnetic stress associated with MHD tur- 
bulence driven by MRI. Following SPB99, we assume that 
the only non-zero components of T are the azimuthal com- 
ponents, 



pr 



d 



(?)■ 



dr \ r 
psin9 d / v,f, 
r 89 Vsin# 



(4) 
(5) 



This is because the MRI is driven only by the shear associ- 
ated with orbital dynamics. Other components of the stress 
are much smaller than the azimuthal components (Stone 
& Pringle 2001). We adopt the coefficient of shear viscos- 
ity p = up. We assume v oc r 1 ^ 2 , which is the usual "a" 
description. In all models, we set a=0.01 to eliminate the 
discrepancy due to the magnitude of viscosity coefficient. 

We use pseudo-Newtonian potential to mimic the gen- 
eral relativistic effects, ip = —GM/ (r — 2r 9 ), where M is the 
center black hole mass, G is the gravitational constant and 
r g = GM/c 2 is the gravitational radius. The self gravity of 
the gas is neglected. 



2.2 Models 

For our aim, we adopt four kinds of initial conditions. In 
all our models, the inner boundaries locate at r = 2.7r 9 . In 
the angular direction, the boundary conditions are set by 
symmetry at the poles. Except the results in Section 3.4, at 
the radial boundaries, we use outflow boundary conditions. 

In model A, the initial condition of our simulation is an 
equilibrium torus with constant specific angular momentum 
given by (Papaloizou & Pringle 1984) 



p _ (7 - 1)GM 
P 



7R0 



Ro 

r 



Ro 
2 V rsmfl 



1 

Id 



(6) 



Here, Ro is the radius of the center (density maximum) of 
the torus, and d is the distortion of the torus. We assume 
Ro = 200r 9 and d = 1.25 in our models. The maximum 
density of the tours p ma x = 1.0. The specific angular mo- 
mentum of the initial tours equals to the Keplerian angular 
momentum at the torus center. Initially, the torus is embed- 
ded in a low-density medium. The density and pressure of 
the medium are p m = 10 -4 and p m — p m /r, respectively. 

In model B, we use the one-dimensional steady global 
ADAF solution (G-ADAF-S; Yuan 1999) as our initial con- 
dition. Because this solution is one-dimensional, we extend 
it into two-dimensions in the following way. The values 
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Figure 1. One-dimensional global solution of hot accretion flows 
on which the initial condition of Model B is based. The dotted 
line in the angular momentum plot shows the Keplerian one; the 
dotted line in the radial velocity plot shows the sound speed. 



of p(r, 0), p(r, 0), v r (r, 0), l(r, 0) are from the global so- 
lution of ADAF. The density and pressure of flow expo- 
nentially decreases along 6: p(r,6) = p(r, 0)exp(— z 2 /2H 2 ), 
p(r,9) — p(r, 0)exp(— z 2 /2H 2 ), where z = r|cos(?|. For sim- 
plicity, we set the radial velocity v r (r,9) = v r (r, 0) and az- 
imuthal velocity v${r,ff) — l(r,0)/(rsm(8)). Initially, the 9 
component of the velocity is set to 0. In Figure 1, we plot 
the radial structure of these initial quantities. 

Model C series are "injection" type. The gas is con- 
tinually injected into the computational domain from the 
outer boundary. The 9 component of the velocity of the 
injected gas is set to be zero. The density of the injected 
gas has a Gaussian distribution about the equator p(r, 9) = 
Pma:rexp[-0.5(6> — 7r/2) 2 ], with p m ax = 1. In model C series, 
we have assumed that the injected gas have different angular 
moment, radial velocity and internal energy. The rotational 
and radial velocities of the injected gas are i;^ = fdiVk and 
v r = f r Vk, respectively, with Vk being the Keplerian velocity 
at the injection radius. The internal energy of the injected 
fluid e = fepip- In models C1-C4, f e is set to 0.2, so the gas 
in these models is marginally bound, i.e. the Bernoulli pa- 
rameter Be < 0. The detailed parameters are listed in Table 
1. 

In model D, the initial condition is adopted to be more 
"realistic". We inject gas in full-range of 9. The properties of 
the injected gas are taken from a parsec-scale, realistic SPH 
simulation of the gas dynamics around Sgr A* (see Cuadra 
et al 2006, 2008 for details). The gas in this simulation orig- 
inates as stellar winds of young massive stars, whose orbital 
and mass-loss properties are well determined from obser- 
vations (Paumard et al 2006; Martins et al 2007). Shocks 
between different stellar winds heat up the gas, creating a 
hot accretion flow. They ran the simulation for ~ 1000 yr, 
which corresponds to several orbital times for the stars that 
originate the gas. In order to apply the 3 dimensional data 
from the simulation by Cuadra (2008) to the current 2 di- 
mensional simulation, we have to chose a symmetry axis - 
we take the vector normal to the plane on which most stellar 
orbits lie as the rotation axis. We then average the data over 
the azimuthal angle to obtain 2 dimension data. Then, the 
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Table 1. Models in this paper 



Run 


Descriptions 


N r X Ng 


l \j<I}) 




e( f 3 ) 

C VJ e J 






s 4 


V 5 


A 


torus 


168 x 88 


0.9-1.1 


0.0 




4.5 


800r 9 


0.49 


0.53 


13 


G-ADAF-S 


200 x 88 


0.4-0.8 


0.01-1 


0.2 


10. 


600r 9 


0.55 


0.8 


CI 


injection 


168 x 88 


0.95 


0.1 


0.2 


900. 


600r 9 


0.53 


0.51 


02 


injection 


168 x 88 


0.55 


0.1 


0.2 


900. 


600r 9 


0.54 


0.49 


C3 


injection 


168 x 88 


0.55 


0.01 


0.2 


900. 


600r 9 


0.55 


0.48 


C4 


injection 


168 x 88 


0.25 


0.1 


0.2 


900. 


600r 9 


0.54 


0.63 


C5 


injection 


168 x 88 


0.1 


0.1 


2.0 


900. 


600r 9 





1.31 


C6 


injection 


168 x 88 


0.1 


0.1 


1.0 


900. 


600r 9 





1.3 


Dt 


injection 


168 x 88 


0.6 


0.01-0.1 


0.1 


600. 


600r 9 


0.47 


0.65 



1,/^ is the ratio of angular momentum to the local Keplerian angular momentum for the initial tours or injected gas. 
2,/r is the ratio of radial velocity to the local Keplerian velocity for the initial tours or injected gas. 
3,/ e is the ratio of the internal energy to the local gravitational energy for the initial torus or the injected gas. 
4, s is the power law index of the radial profile of inflow rate. 
5, p is the power law index of the radial profile of density, 
t The properties of the injected gas in model D are more "realistic". 



2 dimensional data was used as initial conditions in model 
D presented in this work. Looking at these data, we find the 
average specific angular momentum of the gas at the inner 
boundary of their simulation (5000r 9 ) is / ~ 0.6/fc. The ra- 
dial velocity has large fluctuation, and is ~ 0.01 — O.lvk- The 
temperature is ~ 7 x 10 7 K. We use method of least squares 
to fit the results from Cuadra (2008) to obtain the physical 
quantities used at the outer boundary (600r 9 ) in model D 
in the present paper. 

Table 1 summarizes the properties of the models 
mentioned-above. Column (2) gives a brief description of 
the models (initial state is torus, global ADAF solution or 
injection from the outer boundary). Column (3) shows the 
numerical resolution. Columns (4)- (6) show the properties 
of the gas in the initial conditons. They are the specific an- 
gular momentum compared to the local Keplerian angular 
momentum, the radial velocity compared to the local Kep- 
lerian velocity, the specific internal energy compared to the 
local gravitational energy. Column (7) gives the final time 
at which we stopped the simulation (in units of the Keple- 
rian orbital time at r = 200 r g ). Column (8) denotes the 
outer boundary. When the simulations achieve steady state, 
the power- law index of the radial profiles of mass inflow rate 
and density are given in columns (9) and (10), respectively. 



2.3 Numerical method 

We use the ZEUS-2D code (Stone & Norman 1992a,1992b) 
to solve equations (l)-(3). We adopt non-uniform grid in the 
radial direction (Ar)i+i/(Ar)i = 1.037. Similarly, we adopt 
non-uniform angular zones with (A0)j+i/(A9)j — 0.9826 
for < 9 < tt/2 and (A0) j+1 /(A0) j = 1.0177 for tt/2 sC 6 < 



3 RESULTS 
3.1 Torus model 

Following SPB99, we define the mass inflow and outflow 
rates, Mi„ and M out , as follows, 



(GM/C 



Figure 2. The radial profile of the time-averaged and angle in- 
tegrated mass inflow rate M; n (solid line), outflow rate Mout 
(dashed line), and the net rate M acc (dotted line) in model A. 
They are defined in equations (7), JSJ, and J9}, respectively. 



(7) 
(8) 



M in (r) = 27T7- / pmin(u r ,0) sin#d6> 
Jo 

2 H 

M out (r) = 2nr / pmax(i> r , 0) sin8d6 
Jo 

The net mass accretion rate is, 

Macc(r-) = Af in (r) + Mout(r) 

In Figure 2, we plot the time-averaged (from 4-4.5 or- 
bits) and angle-integrated mass accretion rate in model A. 
The radial profile of the inflow rate from 10r 9 to 200r 9 can 
be described by 

/ \ 0.49 

M in = M in (r ut) I — J , (10) 



(0) 



and the description of viscosity in model A are same as those 
in "Run K" of SPB99. The radial profile obtained in SPB99 
is Min oc r 3//4 , steeper than ours. The discrepancy is because 
the coefficient of viscosity in our model is 10 times larger 
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Figure 3. Radial profiles of some quantities in model A. All quan- 
tities have been averaged over time and the polar angle between 
9 = 84° to 6 = 96° . The dotted line in the lower-left plot corre- 
sponds to the Keplerian angular momentum. The dotted line in 
the lower-right plot corresponds to the sound speed. 



than that in "Run K" of SPB99 (IA99; SPB99; Yuan, Wu 
& Bu 2012). 

Figure 3 shows the radial structure of the time-averaged 
(from 4-4.5 orbits) flow near the equatorial plane in model 
A. The solutions are averaged over angle between 9 = 84° to 
8 — 96° . From the inner boundary to 200r 3 , the density, gas 
pressure, angular momentum can be described by a power 
law scaling with radius, 



Figure 4. Same with Fig. 2, but for model B. 



p oc r 



(11) 



This density profile is significantly flatter than that obtained 
in Yuan, Wu & Bu (2012) where they got p oc r -0 ' 85 . This 
is because in that work the initial torus is put at a much 
larger radius, Ro = 10 4 r 9 . Thus, the radial dynamical range 
of Yuan, Wu & Bu (2012) is much larger than the present 
work and the results suffer little from the general relativity 
effect close to the black hole. The density profile is similar to 
that obtained in SPB99, where they found p oc r -0,5 . The 
numerical setting in SPB99 is different from that in this 
paper, SPB99 adopt Newtonian potential and a = 10~ 3 . 
The former inclines to make the density profile steeper while 
the latter inclines to make it flatter. The rough consistency 
between SPB99 and the present work is because the two 
effects cancel each other. 



3.2 Model B 

In model B, our initial condition is based on the one- 
dimensional global ADAF solution (G-ADAF-S). The ra- 
dial structure of our initial conditions is shown in Fig. 1. 
One feature is that the angular momentum of the flow is 
smaller than that of the torus model. Fig. 4 shows the time- 
averaged (from 800-900 orbits) and angle-integrated mass 
accretion rate of model B. The radial profile of the inflow 
rate from 10r 9 to 200r s can be described by 



M in = M in (rw) 



T 'out 



(12) 



(gm/c 2 ; 



(GM/C 2 ) 



(GM/C 2 ) 



(GM/C ) 



Figure 5. Same with Fig. 3, but for model B. 



The result is approximately consistent with Model A. 

Figure 5 displays the radial structure of the time- 
averaged (from 800-900 orbits) flow near the equatorial 
plane (between 9 = 84° to 6 — 96°). From the inner bound- 
ary to 200r g , the density, gas pressure, angular momentum 
can be described by, 



-0.8 

p oc r , p oc r 



I oc r u 



(13) 



We see that the density profile is much steeper than that in 
model A. We can understand it as follows. Because of the 
initial condition, in the outer region (close to 200r 3 ), the 
angular momentum in model A is much higher than that of 
model B; therefore, the radial velocity in model A is much 
smaller than that of model B. Close to the black hole, in 
both models, the radial velocity is equal to light speed. The 
radial profile of velocity in model A is steeper than that of 
model B. Given that mass accretion rate have comparable 
power-law index, the steeper radial velocity profile in model 
A results in flatter density profile. 
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Figure 6. Same with Fig. 2, but for model CI (black line), C2 
(red line), C3 (blue line), and C4 (green line). 



Figure 7. Same with Fig. 3, but for model CI (black line), C2 
(red line), C3 (blue line) and C4 (green line). 



3.3 Injection models 



3.3.1 Injection models with large and moderate angular 
momentum 



We first define the circularization radius r c as 



r c 

Tout 



i 

' o k 



(14) 



where r out is our outer boundary, U n j and Z k are the an- 
gular momentum of the injected gas and the Keplerian an- 
gular momentum at r out , respectively. When l in j <C l ok , the 
injected gas will quickly infall until it get close to the circu- 
larization radius. Because the gas pressure gradient force is 
important in hot accretion flow, the effective circularization 
radius (hereafter, ECR), where the gas stops quick infall, 
may be 2 ~ 3 times bigger than r c . 

Fig. 6 shows the mass accretion rates in models CI ~ 
C4 while Fig. 7 shows the time averaged physical quantities 
near the equator. In model C4, the circularization radius 
r c ~ 37.5r g , Because of the pressure gradient force, the ECR 
of gas locates at ~ 100r 9 . Beyond ECR, we see that the 
angular momentum and inflow rate are both approximately 
constant. The reason why the inflow rate is constant is that 
convection is weak and the accretion timescale is short (refer 
to §3.3.2 for more details). In model C4, quasi-steady state 



50r 3 can be described by 



Min = Min(rw) 



Tout 



(15) 



From the inner boundary to 50r ff , the density, gas pressure, 
angular momentum can be described by, 

-0.63 —1.31 ; 0.2 / 1c n 

peer ,pocr ,( a r . (16) 

The circularization radii of the gas in models CI and C2 
are 541. 5r g and 181. 5r g , respectively. The mass inflow rates 
of models CI and C2 from 10r s to 200r s can be described 



M in cxr a53 ,M in «r - 54 , 



(17) 



respectively. The density, pressure, and angular momentum 
of Model CI and C2 can be described by, 



r(GM/c 2 ) 



Figure 8. Same with Fig. 2, but for model D. 



-0.51 -1.32 , 0.31 

p oc r ,p oc r ,( a r , 

and 

-0.49 -1.31 , 0.33 

par ,p oc r , t oc r , 



(18) 



(19) 



respectively. 

In order to investigate the effects of injected radial ve- 
locity, we carry out model C3. The only difference between 
models C2 and C3 is that the injected radial velocity in C3 
is 10 times smaller than that in model C2. Compared to 
model C2, the mass accretion rate close to the black hole in 
C3 is slightly smaller (by a factor smaller than 2), the profile 
of inflow rate is slightly steeper, while the slops of density, 
pressure and angular momentum in model C3 almost remain 
same. We thus conclude that the effects of radial velocity of 
the injected gas are small. 

In models CI ~ C4, the properties of the injected gas 
are "artificially" given. In Model D , we adopt a "realistic" 
initial condition, i.e., the properties of the injected gas are 
taken from the large-scale simulations of Sgr A* (Cuadra 
2008). Figs. 8 & 9 show the profiles of accretion rates and 
physical quantities. The results are 



M oc r 



(20) 
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(GM/c 2 ) 



Figure 10. Same with Fig. 2, but for model C5 (black line) and 
C6 (red line). 



-0.65 -1.57 , 0.34 

p oc r ,p ocr ,1 ocr 



(21) 



3.3.2 Injection models with very low angular momentum 

Previous analytical works found that when the angular mo- 
mentum of the flow is smaller than a critical value, the ac- 
cretion pattern is Bondi-like, i.e., the sonic radius is large, 
no matter the flow is inviscid (Abramowicz & Zurek 1981) 
or viscous (Yuan 1999). Proga & Begelman (2003, hereafter, 
PB2003) studied inviscid accretion flow and found that when 
the angular momentum of the flow is smaller than the Kep- 
lerian angular momentum at 16r 9 , the gas inside the Bondi 
radius has a constant accretion rate. Because viscosity is 
neglected in PB2003, there is no heating for the gas, the 
entropy of the gas is constant with radius. Therefore the 
flow in PB2003 is convectively stable. When viscous heating 
is included, the entropy of gas should increase inwards and 
convection may play some role on the mass accretion rate 
profile. It is interesting then to investigate whether in this 
case the mass accretion rate is constant. 

In models C5 and C6, the injected angular momen- 
tum of the gas equals the Keplerian angular momentun at 
6r s . In models C5 and C6, the temperature of the injected 



gas is higher than Virial temperature, so the Bondi radius 
Rb = GM/coo (where c x is the sound speed at infinity) is 
smaller than the outer boundary, Rb ~ 300r 9 and 150r 9 , 
respectively. In reality, Rb is usually much larger. 

Figures 10 and 11 show the radial profiles of mass ac- 
cretion rates and some quantities near the equator in models 
C5 and C6. The results are 



M ~ const. 



(22) 



for both models, and correspondingly, the density profile is 
close to the analytical Bondi result, 



-1.31 

p oc r , p oc r 



(23) 



for models C5 and C6. The decrease of mass accretion rate 
close to the outer boundary is not physical but due to bound- 
ary effects. 

We analyze that the reason for a constant accretion rate 
is because convection is very weak. The condition for con- 
vective instability in a rotating accretion flow is: 



N: 



N 



+ a 2 _ IdP d\n(P 



l/-y 



Ip) 



p dR 



dR 



a 2 < 0, 



(24) 



where N is the usual Brunt- Vaisala frequency and a is the 
epicyclic frequency which is equal to rotation angular veloc- 
ity for nearly Keplerian rotation (Narayan & Yi 1994). For 
a non-rotating flow, a = 0, this condition is then equivalent 
to an inward increase of entropy, which is the well-known 
Schwarzschild criterion. In the accretion process, the inward 
increase of entropy is because of viscous heating. The vis- 
cous heating rate is proportional to square of gradient of 
angular velocity. In models C5 and C6, the gradient of an- 
gular velocity is very small. Therefore, viscous heating rate 
is very small in models C5 and C6, which results in a small 
entropy gradient in models C5 and C6. The first term in 
equation (|24|) is proportional to gradient of entropy. We find 
in almost half computational domain the flow is convectively 
stable according to equation (|24[) . In addition, the accretion 
timescale of the flow in models C5 and C6 is very small. 
We find that only in a very small region around 100r 9 , the 
growth timescale of convection is slightly smaller than the 
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gas accretion timescale. This is why convection is unimpor- 
tantQ. 

3.4 Effects of boundary conditions 

In numerical simulations, for the radial boundaries, the out- 
flow boundary conditions are usually adopted. This condi- 
tion means that the hydrodynamical variables (e.g. density, 
velocity, pressure) are directly copied from the first (or last) 
active zone to the ghost zone. In spherical coordinate, the 
mass flux is expressed as r 2 pv r - If we use outflow boundary 
conditions, the gradient of mass flux is not zero at radial 
boundaries, thus there may be artificial mass accumulation 
at the radial boundaries. In order to avoid this problem, 
some authors assume that the gradient of mass flux cross the 
boundaries is zero. We have studied how these two bound- 
ary conditions affect the properties of accretion flow in torus 
model. For outer radial boundary, we find that if the outer 
boundary is much larger than the torus center, the proper- 
ties of the flow are very insensitive to the outer boundary 
conditions. We also find that the properties of the flows are 
also very insensitive to the inner boundary conditions. The 
reason is that at the inner boundary, the flow is supersonic, 
so no wave signal can transfer back to affect the properties 
of the flow. 



4 SUMMARY 

In this paper, we have investigated the effects of initial 
and boundary conditions on the properties of hot accre- 
tion flow using two-dimensional hydrodynamic simulations. 
We have considered several initial conditions, including the 
most widely adopted "torus" and "injection", and the more 
realistic flow properties from large-scale simulation. Special 
attention is paid to the radial profiles of mass accretion rate 
and density. Both profiles can be described by a power-law 
function (M- m oc r s ,/> « r~ p ). The simple description to 
different initial conditions and the corresponding power-law 
index are summarized in Table 1. 

The most interesting result is that for various initial 
conditions, if the angular momentum of the accretion flow 
is not too low so that a rotationally supported accretion flow 
can form beyond 10r 9 , the power-law index of the radial pro- 
file of inflow rate s lies in a narrow range of 0.47 < s < 0.55. 
But on the other hand, the radial profile of density is more 
sensitive to the initial condition. The power-law index p lies 
in a wider range 0.48 < p < 0.8. Given that the inflow rate 
profile is almost same for various initial conditions, the di- 
versity of the density profile is because different initial con- 
ditions give different angular momentum profiles and thus 
different radial velocity profiles. Then it is an interesting 
question why different initial conditions give roughly same 
inflow rate profile. The readers are referred to Begelman 
(2012) for an explanation. 

Since M oc r 2 v r p(r), for hot accretion flows, we usually 

1 In the MHD simulation of Bondi accretion flow, under some ini- 
tial configuration of the magnetic field, Igumcnshchev & Narayan 
(2002) find significant magnetic reconncction heating thus en- 
tropy production. In this case, the flow is convectively unstable 
and p oc r -0,5 . 



expect p = 1.5 — s, if v r oc r -0,5 according to the self-similar 
scaling law of hot accretion flows (e.g. Narayan & Yi 1994). 
We note that our results do not satisfy p — 1.5 — s, but 
smaller than this value. This is partly because the density 
profile depends not only on the inflow rate profile, but also 
on outflow rate profile, which is much steeper. Another rea- 
son is that v r oc r~ ' 5 is not well satisfied since the flow is 
close to the black hole. 

If the angular momentum of the flow is very low (i.e., 
models C5 and C6), we find that inside the Bondi radius, 
the mass accretion rate is almost constant. Physically, this 
is because convection is very weak, thus, very few outflow is 
produced. In these models, the viscous heating rate is very 
small due to the small angular velocity gradient. Small vis- 
cous heating rate results in a very small gradient of entropy, 
thus, convection is weak. In addition, the accretion timescale 
is very small. 

We have also studied the effects of boundary conditions. 
We have considered two widely adopted inner and outer 
boundary conditions, which are "outflow" boundary condi- 
tion and "mass flux conservation" boundary condition. We 
have found that the properties of the flow are very insensi- 
tive to the radial boundary conditions. 
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